A computational study to determine the role of σ-hole in Br/OH substituted nido-carborane and its binding capabilities.

A detailed investigation of the σ-hole on the halogen atom present in the nido-heteroboranes is made by employing quantum mechanical methods. The bromide and the hydroxyl groups are incorporated in the exo-substituents of the nido-boranes. The potential of the bromide σ-hole was compared to that of electrostatic potential of hydroxyl group counterpart. The presence of a carbon atom vertex, in a different position of a system, influences the σ-hole and hence its binding abilities. Bromide substituted nido-carboranes have less potential and hence weaker binding ability compared to their closo-counterparts. Binding affinity with aliphatic is found to be more compared to that of aromatic system. The presence of solvent dampened the electrostatic interactions. Apart from the neutral system, the binding capabilities of charged nido-heteroboranes were also studied. The results of this study will be further useful for several applications viz., crystal engineering, drug designing (Pharmaceuticals), medicine, material science, energy storage devices, etc.

The first HyDRA challenge for computational vibrational spectroscopy.

Vibrational spectroscopy in supersonic jet expansions is a powerful tool to assess molecular aggregates in close to ideal conditions for the benchmarking of quantum chemical approaches. The low temperatures achieved as well as the absence of environment effects allow for a direct comparison between computed and experimental spectra. This provides potential benchmarking data which can be revisited to hone different computational techniques, and it allows for the critical analysis of procedures under the setting of a blind challenge. In the latter case, the final result is unknown to modellers, providing an unbiased testing opportunity for quantum chemical models. In this work, we present the spectroscopic and computational results for the first HyDRA blind challenge. The latter deals with the prediction of water donor stretching vibrations in monohydrates of organic molecules. This edition features a test set of 10 systems. Experimental water donor OH vibrational wavenumbers for the vacuum-isolated monohydrates of formaldehyde, tetrahydrofuran, pyridine, tetrahydrothiophene, trifluoroethanol, methyl lactate, dimethylimidazolidinone, cyclooctanone, trifluoroacetophenone and 1-phenylcyclohexane-cis-1,2-diol are provided. The results of the challenge show promising predictive properties in both purely quantum mechanical approaches as well as regression and other machine learning strategies.

Ranking Power of the SQM/COSMO Scoring Function on Carbonic Anhydrase†II-Inhibitor Complexes.

Accurate prediction of protein-ligand binding affinities is essential for hit-to-lead optimization and virtual screening. The reliability of scoring functions can be improved by including quantum effects. Here, we demonstrate the ranking power of the semiempirical quantum mechanics (SQM)/implicit solvent (COSMO) scoring function by using a challenging set of 10 inhibitors binding to carbonic anhydrase II through Zn2+ in the active site. This new dataset consists of the high-resolution (1.1-1.4 Å) crystal structures and experimentally determined inhibitory constant (Ki ) values. It allows for evaluation of the common approximations, such as representing the solvent implicitly or by using a single target conformation combined with a set of ligand docking poses. SQM/COSMO attained a good correlation of R2 of 0.56-0.77 with the experimental inhibitory activities, benefiting from careful handling of both noncovalent interactions (e.g. charge transfer) and solvation. This proof-of-concept study of SQM/COSMO ranking for metalloprotein-ligand systems demonstrates its potential for hit-to-lead applications.

Synthesis of α-Branched Acyclic Nucleoside Phosphonates as Potential Inhibitors of Bacterial Adenylate Cyclases.

Inhibition of Bordetella pertussis adenylate cyclase toxin (ACT) and Bacillus anthracis edema factor (EF), key virulence factors with adenylate cyclase activity, represents a potential method for treating or preventing toxemia related to whooping cough and anthrax, respectively. Novel α-branched acyclic nucleoside phosphonates (ANPs) having a hemiaminal ether moiety were synthesized as potential inhibitors of bacterial adenylate cyclases. ANPs prepared as bisamidates were not cytotoxic, but did not exhibit any profound activity (IC50 >10 µm) toward ACT in J774A.1 macrophages. The apparent lack of activity of the bisamidates is speculated to be due to the inefficient formation of the biologically active species (ANPpp) in the cells. Conversely, two 5-haloanthraniloyl-substituted ANPs in the form of diphosphates were shown to be potent ACT and EF inhibitors with IC50 values ranging from 55 to 362 nm.

Probing the communication of deoxythymidine triphosphate in HIV-1 reverse transcriptase by communication maps and interaction energy studies.

We calculate communication maps for HIV-1 Reverse Transcriptase (RT) to elucidate energy transfer pathways between deoxythymidine triphosphate (dTTP) and other parts of the protein. This approach locates energy transport channels from the dTTP to remote regions of the protein via residues and water molecules. We examine the water dynamics near the catalytic site of HIV-1 RT by molecular dynamics (MD) simulations. We find that, within the catalytic site, the relaxation of water molecules is similar to that of the hydration water molecules present in other proteins and the relaxation time scale is fast enough to transport energy and helps in communication between dTTP and other residues in the system. To quantify energy transfer, we also calculate the interaction energies of dTTP, 2Mg2+, doxy-guanosine nucleotide (DG22) with their surrounding residues by using the B3LYP-D3 method. The results, from classical vibrational energy diffusivity and QM interaction energy, are complementary to identify the important residues involved in the process of polymerization. The positive and negative interactions by dTTP with different types of residues in the catalytic region make the residues transfer energy through vibrational communication.

Computational study to understand the energy transfer pathways within amicyanin.

Vibrational energy diffusivities between the residues present in Amicyanin copper protein are calculated and presented in form of communication map. From those results energy flow pathways from the copper metal ion to the inter protein residue Glu31 are identified. Our finding suggests many different pathways are possible and copper metal ion in oxidized and reduced state switches the pathways. Our finding also suggests the cooperative nature of surrounding residues and water molecules towards selecting the pathways. The major transport channels in the oxidised state are, Cu2+---> MET28---> LYS29---> TYR30---> GLU31 and Cu2+---> MET98---> TYR30--- > GLU31. And in the reduced state Cu+---> CYS9---> TYR30---> GLU31 and Cu+---> MET28---> LYS2---> TYR30---> GLU31. We studied further the interaction energies between the copper ion and neighbouring residues using B3LYP/QZVP method. Both the methods complement each other in predicting the energy flow pathways and the cooperative nature of residues.

Representative Amino Acid Side-Chain Interactions in Protein-DNA Complexes: A Comparison of Highly Accurate Correlated Ab Initio Quantum Mechanical Calculations and Efficient Approaches for Applications to Large Systems.

Representative pairs of amino acid side chains and nucleic acid bases extracted from available high-quality structures of protein-DNA complexes were analyzed using a range of computational methods. CCSD(T)/CBS interaction energies were calculated for the chosen 272 pairs. These reference interaction energies were used to test the MP2.5/CBS, MP2.X/CBS, MP2-F12, DFT-D3, PM6, and Amber force field methods. Method MP2.5 provided excellent agreement with reference data (root-mean-square error (RMSE) of 0.11 kcal/mol), which is more than 1 order of magnitude faster than the CCSD(T) method. When MP2-F12 and MP2.5 were combined, the results were within reasonable accuracy (0.20 kcal/mol), with a computational savings of almost 2 orders of magnitude. Therefore, this method is a promising tool for accurate calculations of interaction energies in protein-DNA motifs of up to ∼100 atoms, for which CCSD(T)/CBS benchmark calculations are not feasible. B3-LYP-D3 calculated with def2-TZVPP and def2-QZVP basis sets yielded sufficiently good results with a reasonably small RMSE. This method provided better results for neutral systems, whereas positively charged species exhibited the worst agreement with the benchmark data. The Amber force field yielded unbalanced results-performing well for systems containing nonpolar amino acids but severely underestimating interaction energies for charged complexes. The semiempirical PM6 method with corrections for hydrogen bonding and dispersion energy (PM6-D3H4) exhibited considerably smaller error than the Amber force field, which makes it an effective tool for modeling extended protein-ligand complexes (of up to 10,000 atoms).

A comparison of ab initio quantum-mechanical and experimental D0 binding energies of eleven H-bonded and eleven dispersion-bound complexes.

Dissociation energies (D0) of 11 H-bonded and 11 dispersion-bound complexes were calculated as the sum of interaction energies and the change of zero-point vibrational energies (ΔZPVE). The structures of H-bonded complexes were optimized at the RI-MP2/cc-pVTZ level, at which deformation and harmonic ΔZPVE energies were also calculated. The structures of dispersion-bound complexes were optimized at the DFT-D3 level, and harmonic ΔZPVE energies were determined at the same level as well. For comparison, CCSD(T)/CBS D0 energies were also evaluated for both types of complexes. The CCSD(T)/CBS interaction energy was constructed as the sum of MP2/CBS interaction energy, extrapolated from aug-cc-pVTZ and aug-cc-pVQZ basis sets, and ΔCCSD(T) correction, determined with the aug-cc-pVDZ basis set. The ΔZPVE energies were determined for all complexes at the harmonic level and for selected complexes, these energies were also calculated using second-order vibration perturbation (VPT2) theory. For H-bonded complexes, the harmonic CCSD(T)/CBS D0 energies were in better agreement with the experimental values (with a mean relative error (MRE) of 6.2%) than the RI-MP2/cc-pVTZ D0 (a MRE of 12.3%). The same trend was found for dispersion-bound complexes (6.2% (MRE) at CCSD(T)/CBS and 7.7% (MRE) at the DFT-D3 level). When the anharmonic ΔZPVE term was included instead of harmonic one, the agreement between theoretical and experimental D0 deteriorated for H-bonded as well as dispersion-bound complexes. Finally, the applicability of "diagonal approximation" for determining the anharmonic ΔZPVE was shown. For the phenolH2O complex, the ΔZPVE energy calculated at the VPT2 level and on the basis of "diagonal approximation" differed by less than 0.1 kcal mol(-1).

Communication maps computed for homodimeric hemoglobin: computational study of water-mediated energy transport in proteins.

Frequency-resolved communication maps provide a coarse-grained picture of energy transport in nanoscale systems. We calculate communication maps for homodimeric hemoglobin from Scapharca inaequivalvis and sample them to elucidate energy transfer pathways between the binding sites and other parts of the protein with focus on the role of the cluster of water molecules at the interface between the globules. We complement analysis of communication maps with molecular simulations of energy flow. Both approaches reveal that excess energy in one heme flows mainly to regions of the interface where early hydrogen bond rearrangements occur in the allosteric transition. In particular, energy is carried disproportionately by the water molecules, consistent with the larger thermal conductivity of water compared to proteins.

Dynamics of water clusters confined in proteins: a molecular dynamics simulation study of interfacial waters in a dimeric hemoglobin.

Water confined in proteins exhibits dynamics distinct from the dynamics of water in the bulk or near the surface of a biomolecule. We examine the water dynamics at the interface of the two globules of the homodimeric hemoglobin from Scapharca inaequivalvis (HbI) by molecular dynamics (MD) simulations, with focus on water-protein hydrogen bond lifetimes and rotational anisotropy of the interfacial waters. We find that relaxation of the waters at the interface of both deoxy- and oxy-HbI, which contain a cluster of 17 and 11 interfacial waters, respectively, is well described by stretched exponentials with exponents from 0.1 to 0.6 and relaxation times of tens to thousands of picoseconds. The interfacial water molecules of oxy-HbI exhibit slower rotational relaxation and hydrogen bond rearrangement than those of deoxy-HbI, consistent with an allosteric transition from unliganded to liganded conformers involving the expulsion of several water molecules from the interface. Though the interfacial waters are translationally and rotationally static on the picosecond time scale, they contribute to fast communication between the globules via vibrations. We find that the interfacial waters enhance vibrational energy transport across the interface by ≈10%.

Rattling in the cage: ions as probes of sub-picosecond water network dynamics.

We present terahertz (THz) measurements of salt solutions that shed new light on the controversy over whether salts act as kosmotropes (structure makers) or chaotropes (structure breakers), which enhance or reduce the solvent order, respectively. We have carried out precise measurements of the concentration-dependent THz absorption coefficient of 15 solvated alkali halide salts around 85 cm(-1) (2.5 THz). In addition, we recorded overview spectra between 30 and 300 cm(-1) using a THz Fourier transform spectrometer for six alkali halides. For all solutions we found a linear increase of THz absorption compared to pure water (THz excess) with increasing solute concentration. These results suggest that the ions may be treated as simple defects in an H-bond network. They therefore cannot be characterized as either kosmotropes or chaotropes. Below 200 cm(-1), the observed THz excess of all salts can be described by a linear superposition of the water absorption and an additional absorption that is attributed to a rattling motion of the ions within the water network. By providing a comprehensive set of data for different salt solutions, we find that the solutions can all be very well described by a model that includes damped harmonic oscillations of the anions and cations within the water network. We find this model predicts the main features of THz spectra for a variety of salt solutions. The assumption of the existence of these ion rattling motions on sub-picosecond time scales is supported by THz Fourier transform spectroscopy of six alkali halides. Above 200 cm(-1) the excess is interpreted in terms of a change in the wing of the water network librational mode. Accompanying molecular dynamics simulations using the TIP3P water model support our conclusion and show that the fast sub-picosecond motions of the ions and their surroundings are almost decoupled. These findings provide a complete description of the solute-induced changes in the THz solvation dynamics for the investigated salts. Our results show that THz spectroscopy is a powerful experimental tool to establish a new view on the contributions of anions and cations to the structuring of water.